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INTRODUCTION 

The single muon intensity has been measured in the NUSEX experiment at 
various depth up to 10,000 hg/cm 2 s.r.. 

These intensities have been converted to a sea-level muon energy 
spectrum which is used to derive the primary all-nucleon flux. 

From these data we are able to determine the slopes of the primary 
proton and helium spectra, which are thus used in a model for the 
primary composition producing the observed multiple muon rates. 

DATA REDUCTION AND PRIMARY SPECTRUM 

The NUSEX detector is located in the Mt. Blanc tunnel at a vertical 
depth of about 5000 hg cm" 2 s.r. It consists of a cube of 150 t mass 
and 3.5 side, made of 136 horizontal planes 1 cm thick, interleaved with 
planes of tubes of 1 cm x 1 cm cross section, operating in the limited 
streamer mode. More details are reported elsewhere [1]. 

20429 muon events crossing at least 10 layers were recorded in a zenith 
angle range 0° - 75° during an effective working time of 18,946 hours 
from June 1982 to December 1984. The number of events for different 
multiplicities are reported in table 1. 

The single muon intensity at different depths has been found using the 
procedure of Ref. 2. The general relation 

I(h, 6) = I*» K (h) •G ir,K (h, 6)+I P (h) *G p (h, e) (1) 

y y 

has been fitted to the data. Here the ir and K superscripts refer to 
conventional muons from it and K decay, the p superscript to prompt muons 
from charmed particle decays. 
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The angular enhancement functions have been calculated in Ref. 3 . The 
intensity of prompt muons is found < 4% on the considered range of 
depths, so that the second term in (1) has been neglected. The muon 
intensity is reported in Fig. 1, together with the intensity points 
measured with the spark chamber apparatus located in the Mt. Blanc tunnel, 
garage 27. The agreement is excellent, so giving a unique intensity-depth 
set of measurements from 3900 to 10000 hg/cm 2 s.r. very well represented 
by the relation 

I* ,K (h) = ( 7. 63 ±.48) *10 ‘ 7 exp [-h/ (810.44±. 84) ] cm ' 2 s ‘ ‘sr ‘ 1 
Following the procedure of Ref. 2 we derive the primary all-nucleon 
spectrum in the relevant energy range 10 1 3 t 2 , 10 1 ‘* eV : 


dN 
dE o 


(h Q+ ow-2 . 79±. 03 -2 -1 -1 - 

14. 82:. z;E 0 cm s sr GeV 


The rates of events with exacty n muons have been calculated following 
the procedure described in Ref. 4. 

-Yi(A) 

*n = k K i (A) * E ° P n (E 0 ,A)dE 0 

where (E 0 ,A) is the probability to sample n muons of a shower 
resulting from the interaction of a primary with energy E 0 , and mass 
A, and the primary spectrum is described as a superposition of single 
power spectra 

dN - ... “ Y i (A) . . 

dE 0 = H k i (A) E ° (i=p,a,CN0,Mg,Fe) 

Inputs to (Eq,A) come from the energy and radial muon distributions 
calculated by Gaisser and Stanev (Ref. 5) for our experimental site. We 
assume a value of -2.79 for the spectral index of proton and helium nuclei 
as determined from our data and normalize the flux at 10 TeV to the JACEE 
data. The spectral index for CNO and Mg groups is chosen slighty flatter 
as suggested by direct measurements to which we refer for relative 
normalization. The iron spectrum is normalized to a flux of 2.42 10’ 5 
nuclei/ (m 2 *sr*s •GeV/nucleon) at 100 GeV/nucleon, while its spectral index 
is considered as a free parameter. 
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Beyond a rigidity R c = 2 10 6 GeV/c, all spectra steepen to y = 3.0. 

The breaks for nuclei different from protons are at a total energy 

(A/2)R . 
c 

Fig. 2 shows that the measured rates are well described by an iron 
spectral index in the range 2.6 - 2.7. 

This model describes very well both the all-nucleon flux in the range 
10 1 3 - 10 il * eV and the all-particle flux between 10 1 3 - 10 1 6 eV 
(see for exemple the compilation of Hillas in Ref. 6). Normalization 
coefficients, slopes and breaking points or all components are summarized 
in table 2. 


CONCLUSION 

We have determined a primary cosmic ray spectrum fitting both our 
experimental multiple muon rates and the all-nucleon flux derived from 
the single muon intensities underground. 

In the frame of the interaction model developped by Gaisser, Elbert and 
Stanev, we are able to reproduce NUSEX muon data with a primary 
composition in which the iron spectrum is only slightly flatter than the 
proton one. 

This result rules out the popular idea that the primary composition 
varies drastically with increasing energy, leading to the dominance of 
heavier nuclei at energies 10 1 5 - 10 16 eV. 
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Figure Captions 

Fig. 1 - Muon intensity underground at Mt. Blanc. 

Fig. 2 - Comparison between experimental rates and the predictions of 
our model of primary composition for different values of the 
iron group spectral index. 
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Rate of multiple muons 


Multiplicity 


# events 

Time 



i 


20429 




2 


211 




3 


29 

18946 hours 



4 


4 




5 


2 




6 


1 






Table 2 




Group 

K (m" 2 s ' 

*sr ‘ l GeV ‘ *) y 


Points 

of 


i 

i 

change 

of slope, 

E(GeV) 

P 

3.28 10 1 * 

2.79 

2 

10 6 


He 

1.75 10- 

2.79 

4 

10 6 


CNO 

6.20 10 3 

2.71 

1.4 

10 7 


Mg 

9.20 10 3 

2.71 

2.6 

10 7 


Fe 


2 . 6 t2 . 7 

5.2 

10 7 




h (hg cm 2 ) 




